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Introduction: 

In the late 80 ’s, conceptual designs for an in situ oxy- 
gen production plant were documented in a study by 
Eagle Engineering [1]. In the “Summary of Findings” 
of this study, it is clearly pointed out that: “reported 
process mass and power estimates lack a consistent 
basis to allow comparison.” The study goes on to say: 
“A study to produce a set of process mass, power, and 
volume requirements on a consistent basis is recom- 
mended.” Today, approximately twenty years later, as 
humans plan to return to the moon and venture be- 
yond, the need for flexible up-to-date models of the 
oxygen extraction/production process has become 
even more clear. 

Multiple processes for the production of oxygen from 
lunar regolith are being investigated by NASA, aca- 
demia, and industry. Three processes that have shown 
technical merit are molten regolith electrolysis, hydro- 
gen reduction, and carbothermal reduction. These 
processes have been selected by NASA as the basis for 
the development of the ISRU System Model Tool 
(ISMT). In working to develop up-to-date system 
models for these processes NASA hopes to accomplish 
the following: (1) help in the evaluation process to 
select the most cost-effective and efficient process for 
further prototype development, (2) identify key pa- 
rameters, (3) optimize the excavation and oxygen pro- 
duction processes, and (4) provide estimates on energy 
and power requirements, mass and volume of the sys- 
tem, oxygen production rate, mass of regolith required, 
mass of consumables, and other important parameters. 
Also, as confidence and high fidelity is achieved with 
each component’s model, new techniques and proc- 
esses can be introduced and analyzed at a fraction of 
the cost of traditional hardware development and test 
approaches. A first generation ISRU System Model 
Tool has been used to provide inputs to the Lunar Ar- 
chitecture Team studies. 

Model Description: 

A typical end-to-end ISRU system model is composed 
of a regolith excavation system, regolith feed system, 
chemical processing plant, and liquefaction and stor- 
age system. The system model is divided into modules 
that represent unit operations (e.g., electrolyzer, gas 
separator, reactor, liquefaction, etc). This modularity 
(plug-n-play) feature allows the use of the same unit 
operation model in different oxygen production sys- 
tems simulations, resulting in comparable and consis- 


tent results. Each unit operation is modeled theoreti- 
cally using Excel and Visual Basic for Applications 
(VBA), and validated using available experimental 
data from on-going laboratory work. Each module 
contains a worksheet called “Databus” that functions 
as an interface between modules. This Databus con- 
tains input and output fields, where parameters are 
grouped into arrays using the Range Name feature of 
Excel. These arrays are m- by-1 matrices (m being the 
number of parameters in the array), and can include 
strings, booleans, and numerical values. 

Modules are linked to each other using ‘flow arrays’ 
containing temperature, pressure, and flow rates of 
each compound present in the stream. Furthermore, 
each module contains “Global” Inputs/Outputs and 
“Design” Inputs/Outputs. Global I/O are those pa- 
rameters of interest at the high level such as location of 
lunar outpost, type of power system, mass, power, vol- 
ume, etc. Design I/O are those parameters that are 
specific to each module such as efficiency, diameter, 
materials, etc. 

ISRU System Model Tool: 

The ISMT consists of sub-systems integrated using a 
commercial off-the-shelf model integration software 
[2]. This software offers a graphical interface to link 
or connect each Excel model file, a trade study tool, a 
parametric study tool, and an optimization tool. The 
ISMT consists of: Excavation sub-system, Regolith 
Handling sub-system, Reactor sub-system, Electro- 
lyzer sub-system, Liquefaction sub-system, and Ther- 
mal Energy sub-system. The following is a brief de- 
scription of the capabilities and features of these sub- 
systems: 

Excavation sub-system: 

• Force module calculates forces on digging tools 
and wheels/tracks based on classical soil mechan- 
ics correlations; dimensions of digging tool, 
wheels, and chassis; and power/energy require- 
ments for digging and driving operations. 

• Mass module calculates the mass of individual 
components (digging tool, boom arm, motors & 
actuators, chassis, etc.) based on the dimensions 
and forces calculated in force module. 

• Options for bucket wheel, front-end loader, back- 
hoe, and bull-dozer blade. 

• Options for wheels or tracks. 

• Options for continuous or intermittent digging. 
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• Capability to vary vehicle velocity, depth of cut, 
soil properties, surface slope (for driving), digging 
rate (e.g., regolith per day), per-delivery load, 
down-time between deliveries, and many other 
variables. 

Regolith Handling sub-system: 

• Includes feed and dump hoppers and augers. 

© Hoppers are sized based on amount of regolith 
stored. 

• Augers have an option for heat exchange between 
the feed auger (cold regolith) and dump auger 
(hot-spent regolith). 

Reactor sub-system: 

• Options for carbothermal processing, hydrogen 
reduction processing, and molten regolith elec- 
trolysis. 

• Carbothermal processing option includes a car- 
bothermal reactor, desulfurization unit, and 
methanation reactor. Carbothermal reactor, desul- 
furization, and methanation reactor models are 
based on recent design concepts [3]. 

• Hydrogen reduction processing option includes a 
hydrogen reduction reactor and desulfurization 
unit. Hydrogen reduction reactor has options for a 
fluidized bed, loosely-packed bed, and a rotating 
bed (currently being developed) [4]. 

Electrolyzer sub-system: 

• Options for proton exchange membrane (PEM) 
and solid oxide (SO) electrolyzer. 

• PEM electrolyzer option includes a micro-channel 
heat exchanger, condenser, water pump, phase 
separator, PEM electrolyzer, and a water absorp- 
tion bed. 

• SO electrolyzer option includes gas phase separa- 
tors and SO electrolyzer. 

Liquefaction sub-system: 

• The liquefaction sub-system includes radiators, 
cryocoolers, and storage tanks. This sub-system 
has an option for single or multiple cryocoolers 
for condensing the oxygen produced and manag- 
ing boil-off. 

Thermal Energy sub-system: 

• The thermal energy sub-system includes a rigid 
solar concentrator [3] to provide thermal energy to 
the reactors. An inflatable solar concentrator op- 
tion is currently being developed. 

The ISMT is capable of performing optimization on an 
ISRU system composed of one option of each sub- 
system. Furthermore, the ISMT allows for trade stud- 
ies of each system configuration by replacing any op- 
tion within each sub-system. 


The latest version of the ISRU System Modeling Tool 
will be presented along with results and trade studies 
of various system configurations. 
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ISRU System Model Tool 


Capability: 

- System optimization 

- Allows integration of various modeling a pplica tions: Exce 
MathCad 

- Interchangeability df sysfemjpmponents (“Plug-n-Play”) 

- ISRU systems are easily updated 


cations: Excei^MatLab, and 


Functionality: -nra m 

- Evaluation and of different 02 prodiJrction processes and Excavation 

systems: mass, power, volume V; 

• • 

- Identify key parameters \J' y 

- Evaluate system performance 

- Identify design needs and improvements 


Utility: 

- Prove feasibility of ISRU 

- Support mission architecture studies 

- Support hardware development 






Description of oxygen production processes: 
Carbot he rmal Re duction - PEM Electrolyzer 
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Description of oxygen production processes: 
Ho R eduction - PEM Electrolyzer 
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Description of oxygen production processes 
H 2 Reduction - SO Electrolyzer 
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Description of oxygen production 
processes: Electrolyzer Sub-system 




H20, H2 In Anode-fed PEM 











System Model: A Modular Approach 

ISRU System Level 


Mission Requirements: 

•Number of missions 
•Number of EVAs 

•Description of outpost 


Mission Flequirements 

ISRX 




lnputs: f : 

•Regolith excavation rate 
•Location of excavation site 
•Location of 0 2 plant 
•Type of excavator/transporter 


Inputs: 

•Oxygen production rate 
•Location of 0 2 plant 
•0 2 production process 


0 2 Production 
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Outputs: 

•Regolith required 
•System mass 
•System power 
•System volume 


Outputs: 

•System mass 
•System power 
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System Model: A Modular Approach 

0 2 Production and Excavation System Level 


Oxygen production system is 
divided into modules or unit 
operations 

Excavation system is divided into a 
force module and mass module 

Each module is modeled in a stand-" 
alone Excel workbook using Visual 
Basics for Application (VBA) 

An Excel worksheet functions as the 
Input/Output interface (Databus) 

A VBA module is used to build a 
‘master’ function where the 
calculations are performed 

Modules are linked using Phoenix 
Integration ModelCenter 

Each link represents information 
passed from one module to the 
other 

- The information is passed in the 
form of 




Module 

Molten Regolith 
Electrolysis 

H2 Reduction 

Hopper 

X 

X 

Regolith Handling 

X 

X 


Electrolysis cell X 


.j* •will n. J...*:. _ .jt i i ru.j.i »n f 

Carbothermal 

Components 

X 

Hopper 

X 

Auger, auger-heat 
exchanger 


Reactor chamber, 
electrode, cathode 


Carbothermal ^ Reactor chamber, 

reduction reactor rake system 




• individual cell 

• array of cells 

- Each cell or groups of cells is given 
a specific name (aka: Named ~ 
Range) 


PEM and solid oxide 

Water Electrolyzer X X electrolyzer, water 

pumps 


Water removal unit 


X 

X 

Water removal unit 

De-sulfurization unit 

X 

X 

X 

De-sulfurization unit 


Solar concentrator 


X 


X 


Solar concentrator, 
fiber optics 




System Model: A Modular Approach 

0 2 Production System Level cont’d 


There is a hierarchy in the Databus: 

• Inputs: 

Globallnput: all constants and quantities that are 
specified at the system level 

> Location 

■ > Required 02 Production 

^ ► Designlnput: all other input parameters required to 
run a component but are not generated byootfier 
components 

> Vessel diameter ' 

■ > Material of construction 

^ P Inflow: input values that describe fluid stream 
condition and composition 

> Temperature 

> Pressure 

■ > Flow rates (H 2 , H 2 0, 0 2 , CH 4 , C0 2 , CO, H 2 S, Ar) 

^ P Interface: all other input values a component requires 

that come from another component 

> Reaction Time 

> Batch Time 

> Initial Soil Temperature 

• Outputs: 


GlobalOutput: all calculated values pertinent to 
overall model’s conclusions 

> Mass 

> Power 

> volume 

•Outflow: output values that describe fluid stream 
condition and composition 

> Temperature 

> Pressure 

> Flow rates (H 2 , H 2 0, 0 2 , CH 4 , C0 2 , CO, H 2 S, Ar) 

• DesignOutput: all other calculated output values that 
describe component specifics but are not required by 
other components 

> Vessel height - • 

> Element mass 


r ~ 

B C D E _ ! F 

j 3 

Input Ffld • - ” • - 1 

4 

Value Units 

Description 


5 

Shackleton Rim (short ni< 

ghts) 

Location of Power System 

Modulelnput 

Globallnputs 

6 

Solar (day only] 


Power Type 

7 

1.63000 

m/s2 

Gravity Constant 

8 

Input Solids per Batch 


Flag to specify primary input driver 

3 

15.00000 

*2 

Solid mass in reactor per batch 

10 

1000.00000 

*2 

Desired 02 production per year 

11 

0.80000 


Regolith particle sphericity 

12 

0.00000 

m 

Smallest particle size 

13~ 

000020 

m 

Largest particle size 

14 

3000.00000 

kq/m3 

Average density of solid 

15 

500000 

V. 

FeO percent of solid by mass 

16 

0.20000 

V. 

Sulfur percent of solid by mass 

17 

050000 

hr 

Time to heat up 1 batch of regolith 

18 

12.00000 

batches/day 

Number of batches run per 24 hour day 

19 

Al 2014-T6 


Material of Construction for Reactor and Cyclone 

Designlnputs 

20 

2.00000 


Safety factor for Reactor shell 

l~21 

2.00000 


Safety factor for cyclone 

22 

0.16400 

m 

Diameter (inside) of reactor 

23 

0.02000 1 

m 

Insulation thickness 

24 

Fluidized 


Type of reactor bed 

25 

| 26 

0.30000 


Void fraction in fixed bed 

0.40000 


Void fraction in inciepently fluidized bed 

27 

0.50000 


Void fraction in fluidized or loosely packed bed 

28 

15.00000 

mfs 

Cyclone inlet velocity(15 mfs recommended) 

29 

425.00000 

K 

Initial solid temperature 

Interf aceln 

30 

1100.00000 

K 

Inlet temperature of gases 

Modulelnflov 

31 

206843.00000 

Pa 

Inlet pressure of gases 

32 

0.04800 

mol/s 

h2 molar flowrate in 

33 

0.00000 

mol/s 

h2o molar flowrate in 

34 

0.00000 

mol/s 

o2 molar flowrate in 

35 

36 

0.00000 

molfs 

ch4 molar flowrate in 

0.00000 

mol/s 

co2 molar flowrate in 

37 

0.00000 

mol/s 

co molar flowrate in 

38 

0.00000 

mol/s 

h2s molar flowrate in 

39 

0.00000 

mol/s 

inert molar flowrate in 

40 






41 

Value 


Description 

. wmmnrsmnTrm 

42 

1.15929 

m 

Vessel Height 

ModuleOutput 

DesignOutput 

43 

0.16400 

m 

Vessel Diameter 

44 

009086 

m 

Heater Rod diameter 

45 

1.62182 

kg 

mass of reactor structure 

46 

0.02449 

m3 

volume of reactor structure 

47 

3.57008 

*2 

mass of insulation 

48 

0.01471 

m3 

volume of insulation 

49 

0.02350 

kg 

mass of the cyclone structure 

50 

0.00044 

m3 

volume of the cyclone structure 

51 

0.04670 

mol/s 

required flow rate of hydrogen into reactor 

52 

15.00000 

kg 

required regolith in reactor per batch 

53 

1.18013 

hr 

required time to react each batch of solids 

54 

0.50000 

hr 

required time to heat solids to reaction temperature 

55 

0.87856 

KW 

insulation heat loss 

56 

0.02865 

KJ 

endothermic heat loss 

57 

11562.22982 

KW 

energy required per batch to heat up solids 

58 

1100.00000 

K 

outlet temperature 


ModuleOutflow 

59 

176317.56203 

Pa 

outlet pressure 

60 

0.04427 

molfs 

h2 molar flowrate out 

61 

000221 

mol/s 

h2o molar flowrate out 

62 

0.00000 

mol/s 

o2 molar flowrate out 

63 

0.00000 

mol/s 

ch4 molar flowrate out 

64 

0.00000 

mol/s 

co2 molar flowrate out 

65 

0.00000 

mol/s 

co molar flowrate out 

66 

0.00022 

mol/s 

h2s molar flowrate out 





module 


System Model: A Modular Approach 

0 2 Production Sub-System Level 

The Named Range creates a common interface for 
the modules 
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ISRU System Model Tool: 
Sub-systems Description 

Excavation sub-system: 

- Force module calculates forces on digging tools and wheels/tracks based on classical soil mechanics correlations; 
dimensions of digging tool, wheels, and chassis; and power/energy requirements for digging and driving operations. 

- Mass module calculates the mass of individual components (digging tool, boom arm^otors & actuators, chassis, 
etc.) based on the dimensions and forces calculated in forcejiiodule. 

- Options for bucket wheel, front-end loader, backhoe, and bull-dozer bla^le 

- Options for wheels or tracks. 

- Options for continuous or intermittent di| 
slope (for driving), digging rate (e.g., regSfittfper day), per-delivery load, do 
other variables. 

Regolith Handling sub-system: 

- Includes feed and dump hoppers and augers. 

- Hoppers are sized based on amount of regolith stored* 

- Augers have an option for heat exchange between the fi 
regolith). 

Reactor sub-system: 

- Options for carbothermal processing, hydrogen reduction processing, and molten regolith electrolysis. 

- Carbothermal processing option includes a carbothermal reactor, desulfurization unit, and methanation reactor. 
Carbothermal reactor, desulfurization, and methanation reactor models are based on recent design concepts. 

- Hydrogen reduction processing option includes a hydrogen reduction reactor and desulfurization unit. Hydrogen 
reduction reactor has options for a fluidized bed, loosely-packed bed, and a rotating bed (currently being developed). 

Electrolyzer sub-system: 

- Options for proton exchange membrane (PEM) and solid oxide (SO) electrolyzer. 

- PEM electrolyzer option includes a micro-channel heat exchanger, condenser, water pump, phase separator, PEM 
electrolyzer, and a water absorption bed. 

- SO electrolyzer option includes gas phase separators and SO electrolyzer. 

Liquefaction sub-system: 

- The liquefaction sub-system includes radiators, cryocoolers, and storage tanks. This sub-system has an option for 
single or multiple cryocoolers for condensing the oxygen produced and managing boil-off. 

Thermal Energy sub-system: : 

- The thermal energy sub-system includes a rigid solar concentrator to provide thermal energy to the reactors. An 
inflatable solar concentrator option is currently being developed. 


i auger (cold regolith) and dump auger (hot-spent 
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ISRU System Model Tool: 
Summary of Inputs 



System model optimization: 

- Objective Definition: 

• minimize system mass 

- Design Variables: 

• Mass of regolith per batch 

• Time to heat up or 

• Reactor diameter 

• Number of batches per day 

Location: 

- Equatorial & Polar region 

Excavation sub-system: 

- Single scoop front-end loader, 6 wheels, 'sand combined dig/haul vehicle(s) 

- Intermittent digging (i.e., downtime between deliveries) selected to minimize total 
energy 

Regolith Handling sub-system: 

- No heat exchanger 

Reactor sub-system: 

- Carbothermal processing 

- Hydrogen reduction processing: 2-fluidized bed for continuous operation 

Electrolyzer sub-system: 

- Proton exchange membrane (PEM) 

Liquefaction sub-system: 

- 1 -cryocooler 

- 2-storage tanks 

Thermal Energy sub-system: 

- Rigid solar concentrator 



ISRU System Model Tool 
System Mass Results 
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• All ISRU system mass show a liner dependency with respect to oxygen 
production 

• Carbothermal @ Polar < @ Equatorial: 

- No dependency on regolith composition 

- Operating: Polar = 255 days/yr vs Equatorial = 183 days/yr 

• H2 reduction @ Polar > @ Equatorial: 

- Dependency on regolith composition: FeO =15% equatorial vs 5% polar 

- Even though operating time in equatorial region is less, the effect of regolith 
composition on the overall system mass is greater than the operating time 

• Mass payback is achieved after the first six months of deployment on either 02 
production system 
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ISRU System Model Tool: 
Mass Breakdown - Polar Region 
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□ Regoltih Handling Sub-System 

■ Liquefaction and Storage Sub-system 

□ Thermal Energy Sub-system 

□ Eelctrolysis Sub-System 

■ Reactor Sub-system 

□ Excavation Sub-system 


Numbers on the columns represent the mass of the sub-system in kilograms 


Inefficiency of H2 reduction is seen on the mass of excavation sub-system and 
regolith handling sub-system 

Mass of thermal energy, electrolysis, and liquefaction sub-systems in H2 
reduction is reduced because multiple reactor sub-system is assumed allowing 
continuous processing (no downtime) 

As oxygen production rate increases the fidelity of the reactor sub-system model 
decrease 

- More noticeable in carbothermal model 

- Available experimental data is targeted to low production rate 




ISRU System Model Tool: 

Mass Breakdown - Equatorial Region 



H2 Reduction - Equatorial Region Carbothermal Reduction - Equatorial Region 
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□ Regoltih Handling Sub-System 

■ Liquefaction and Storage Sub-system 

□ Thermal Energy Sub-system 

□ Eelctrolysis Sub-System 

□ Reactor Sub-system 

□ Excavation Sub-system 


Numbers on the columns represent the mass of the sub-system in kilograms 


• Downtime at Equatorial region is greater than at Polar region 

- The effect of operating time is greater in carbothermal than in H2 reduction due to 
multiple reactors on H2 reduction system 




Summary & Conclusions 


ISRU System Model Tool: 

- has been develop for the analysis of ISRU systems: Excavation^nd 02 production 

- is flexible, and allows any component mo^^developed in various applications to be 

integrated into the model environment as long^is the model meets interface 
requirements Ml 


X xr^T N m 

The ISRU System Model Tool hais been successfully used to provide ISRU 
system mass, power, and volum%£pmates to architecture studies 


x . ‘ iX .‘‘m 


Models are validated with available experimental data at low production rates; 
experimental data at higher production rates is needed. 


ISRU System Model Tool mass estimates results show that mass payback is 
achieved within six months of operation 
- Including H2 reduction with Mare-type regolith (-5% FeO) 
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